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We show that the one-particle spectral functions of the one-dimensional Hubbard model diverge at 
the Fermi energy like \uj — ef\~ 3 ^ 8 in the U — > +cxd limit. The Luttinger liquid behaviour \u — ep\ a , 
where a — > 1/8 as U — > +oo, should be limited to \u — ef\ ~ t 2 /[7 (for {/ large but finite), which 
shrinks to a single point, u> = ejr, in that limit. The consequences for the observation of the 
Luttinger liquid behaviour in photoemission and inverse photoemission experiments are discussed. 

(January 18, 1995) 
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Due to large quantum fluctuations, the low-energy 
physics of interacting electrons in ID is not of the Fermi 
liquid type, but can be described by the Luttinger liquid 
theory jl]||]|]. According to that theory, the momen- 
tum distribution function should have no step at kp but 
should behave like n k — n kp cx sign(fc F — k)\k — k F \ a , 
where a is a non-universal exponent that depends on the 
interaction. This has been confirmed for the Hubbard 
model by Ogata and Shiba Q in the U — > +oo limit, 
in which case a — 1/8 ||. Similarly, the local (momen- 
tum averaged) one-particle spectral functions A(u>) (in- 
verse photoemission spectrum) and B(u>) (photoemission 
spectrum) defined by 

A{w) =£|</,iV + l|a+ iCT |0, N) J' 6(u-E? +l +E?) 

f,<r 

B(oj)=^2\(f,N-l\a 0!> 

f,<* 

are expected to behave like \w — efI" close to the Fermi 
energy. In recent photoemission experiments, the func- 
tion B(u) has been measured over a large energy range 
with a resolution of several meV J(|. Under these con- 
ditions, is it possible to detect the Luttinger liquid be- 
haviour? The Luttinger liquid theory itself does not 
predict how far the power law behaviour will hold, and 
to answer this crucial question, a determination of the 
spectral function for microscopic lattice models is nec- 
essary. This is a very difficult problem. Even for the 
Hubbard model, which is soluble by Bethe ansatz (?]] and 
whose low-energy properties are reasonably well under- 
stood , an exact calculation of the spectral functions 
has not been possible so far, although much has been 
done in this direction, both analytically using bosoniza- 
tion JTo| , canonical transformation |n[] and Bethe Ansatz 
|L2j , or numerically with Monte Carlo calculations |l|| . 

In this Letter, we present a calculation of these spectral 
functions for the Hubbard model in the U — > +oo limit. 
The Hubbard model is defined by 



H 



i+l,a a i,a 



h.c 



U 



(2) 



where t is the hopping integral and U is the on site repul- 
sion. We will denote by N the number of fermions and by 
L the number of sites (we choose L even), the site index 
i runs from i — to L—X. Furthermore, we take TV to be 
of the form An + 2 (n integer) so that the ground-state 
is non-degenerate. In the U — > +oo limit, it has been 
shown, using the Bethe Ansatz solution, that the eigen- 
states can be written as a product of a spinless fermion 
wave-function and a squeezed spin wave-function B : 



\NJ) = \i>l Q ({I}))®\x7(QJQ)) 



(3) 



The spin wave-function \x) is characterized by the num- 
ber of down spins N±, the total momentum Q, and 
the quantum number fg within the subspace of mo- 
mentum Q. The spinless fermion part is an eigen- 
function of N spinless fermions on L sites with momenta 



+ Q, where the Ij 



j, j = 1 . . . N, are integer 
quantum numbers. The two components of the wave- 
function are coupled through the momentum Q of the 
spin wave function, which imposes a twisted boundary 
condition on the spinless fermion wave-function (each 
fermion hopping from site L — 1 to site will acquire a 
phase e J< 2). 

In the limit U — > +00, all the states with different spin 
configurations are degenerate and the energy is equal to 
— It coskj, i.e. it does not depend on the quantum 

numbers Jq. In the ground-state, the spinless fermion 
wave-function \^qs l ) is described by the quantum num- 
bers Q = it and {1} = {-N/2, N/2 - 2, N/2 - 1}, so 
that the distribution of the kj 's is symmetric around the 
origin, the spin part being the ground-state of the Heisen- 
berg model according to Ogata and Shiba's prescription 

I- 

For a less than half-filled model, the spectral function 
A(u>) has contributions from both the lower and upper 
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Hubbard bands and A{u) = A lhb (lu) + A UHB (u). 
If we use the Ogata-Shiba wave function into Eq. ([I]), 
we get A lhb (uj) and B(ui). To do this, we just write 
operator aj ^ entering Eq. (|l|) as b\z\ a , where fot creates 
a spinless fermion at site j and zj CT inserts a spin cr 
after skipping the first j spins. A UHB (us) can actually be 
obtained in that framework by a minor modification of 
the wave- function []l4|| . The summation over the different 
spin configurations (/q) can then be performed because 
the energy depends only on the quantum numbers {/} 
and Q, and we get 



A LHB (u;) 
A UHB (u) 



J2C*,n(Q)A q (u;) 

Q.cr 

1— Yl D n ^(Q)B q _q(U- 



L-N- 



>,Q,<T 



J2d^n(Q)B q (u 



(4) 



Cv,n(Q) and D a ^ N (Q) are given by 



Ca l N(Q)=T / \(XN+l(Q,fQ)\Zl <7 \XN(7T,0))\ (5) 

/q 

with Q = 2nj/(N + 1) (j integer) and 

AxJv(Q) = J2 \(XN-i(Q,f Q )\Z MXN(ir,0)}\ 2 (6) 

with Q = 2nj / (N — 1), while Aq(uj) is given by 

E|«^(m)i^i^5,L>| 2 ^-^f +1 +o (7) 
{i} 

and a similar definition holds for Bq(oj). 

The problem has now been reduced to the calculation 
of quantities involving only the charge or the spin. A sim- 
ilar approach has already been followed by Sorella and 
Parola |l5| in their calculation of the momentum distri- 
bution function. For instance, our D a ^{Q) is equivalent 
to their Z(Q). However, they could not calculate the 
spectral function away from half-filling because they did 
not know how to evaluate the charge part in that case. In 
the following, we show how to calculate these quantities 
for any band-filling. 

Let us start with C CTj at((5) and D a ^M{Q)- These 
quantities satisfy the sum rules C<j.n{Q) = 1 and 
J2q D,j,n(Q) = N a /N and they have a singularity at 
Q a = ttNct/N. We have calculated them for small clus- 
ters (Fig. Q), and we found that for = iVj_, in which 
case Q] = Qi = 7r/2 = Qo, they behave like ~ \Q—Qo\~ v 
with X] = 0.49 ± 0.01 if Q > Q and r] = 0.14 ± 0.01 if 
Q < Qo for C Q and 77 = 0.49 ± 0.01 if Q < Q and 
j] = 0.34 ± 0.02 if Q > Qo for D Q . Let us note that, 



although they are very small, C a ^{Q) and -D CT ,./v(Q) do 
not vanish identically for Q < Qo and Q > Qo respec- 
tively, contrary to what was claimed in Ref. ]l5| |. The 
exponents of the main singularities are consistent with 
the theoretical value 1/2. In Fig. [j], the solid lines are 
fits to the numerical results that have been used in the 
following. 

Let us now turn to the charge part. For Q = 7r, the 
only excited states that contribute are those with one 
particle-hole pair, and A^(lo) is just the spinless fermion 
density of states l/ny/it 2 — oj 2 . The problem is not so 
simple for Q 7^ 7r because we have to evaluate matrix ele- 
ments between states with different boundary conditions. 
Let us suppose that k and k' correspond to boundary con- 
ditions with Q and Q' respectively. Then, one can easily 
show that: 



«fc'.a fc 



1 -i(k'-k)/2 i(Q'-Q)/2 sin([Q' - Q]/2) 

L ' sin([k' - k]/2) 



(8) 



Clearly, we are faced with Anderson's orthogonality 
catastrophy Jl6[ and states with many particle — hole ex- 
citations will contribute. More generally, the matrix ele- 
ment (0\a kN ... a k2 a ki aiy k , 2 . . . o v |0) is given by 



-NJ(Q'-Q)N/2 



-<(*j-k,)/a sh jv 



• —l k\ — ki . —1 k\—k 2 

sin 1 2 sin 1 2 



sin 



-1 k' 2 —k\ 



sin 



-1 k' 2 -k 2 



sin 



sin 



Q' -Q 
2 

-1 k[-k N 
2 

-1 k' 2 —k N 



• —1 kj^ — ki • — 1 k' N — k 2 

sin " sin " 



-1 k N — 



The central observation is that this determinant is very 
similar to the Cauchy determinants |]l7| and that it can 
be expressed as a product: 



sin 



j>i 



3>i 



sin - 



n 



Sill 



-1 K 



where the sign is + for N = 1, 4, 5, 8, 9, .. and — for N — 
2, 3, 6, 7, ... After a straightforward calculation one finally 
gets: 



n 



2 _ 1 2N Q 

- L 2N+1 C0S 2 



X 1 1 Sm 

j>i 



2 kj fej j-J 



2 "3 '°i TT • -2 ,XJ 3 



V — V 
j i TT 

sm {I 



sm 



(9) 



For the matrix element entering Bq(uj), a similar ex- 
pression holds with N replaced by N — 1. With such 
an expression, it becomes possible to calculate the spec- 
tral functions numerically. To do that, we generate the 
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quantum numbers Ij and calculate the energy and the 
product of Eq. (||). It turns out that the sum rules 
J™ x A Q {uj)duj = 1-N/L and B Q {uj)duj = N/L hold 
separately for each Q. Moving from Q = it to Q = 0, the 
effect of the many-particle excitations becomes impor- 
tant. The weight of the Van-Hove singularity decreases; 
a power-law singularity appears near the Fermi energy 
and a tail appears beyond the Van-Hove singularity. This 
tail comes from the incoherent part of Aqik ~ tt, uj) and 
Bcj{k ~ 0, uj), and its weight increases when Q decreases. 
Including up to three particle-hole excitations, we al- 
ready get almost all the spectral weight. For instance, for 
L = 60 we get 99.993% of the total weight in the worst 
case, namely for Q = 0. So, for all practical purposes, 
we can limit ourselves to three particle-hole excitations. 
This is an important observation because the number of 
states grows exponentially with the size, and to take into 
account all the states is possible only for smaller systems 
(L ps 30). 

We are now in a position to calculate numerically A(uj) 
and B{ui) and a typical result is presented in Fig. |2|. 
Note that the following sum-rules are satisfied Jl8| , pT| : 
JZo A LHB {Lo)dio = 2(1 - N/L) and A UHB (uo)duo = 
SZ o B{uj)duj = N/L / 

The most interesting and surprising result is that, in- 
stead of going to zero, as the Luttinger liquid theory 
predicts, the spectral functions increase when u — ► Ep, 
where Ef — —2t cos2kF is the Fermi energy of the spin- 
less fermions and hp — ttN/2L is the Fermi momentum 
of the fermions with spin. In fact, using the same frame- 
work, we can prove analytically that they diverge in that 
limit. To do that, we first use the fact that the low en- 
ergy spectrum consists of towers of excitations centered 
at momenta (N + 1)Q/L + Apk F , p = 0,±1,±2, ... to 
write Aq (uj) as the sum of the contributions coming from 
each tower Aq(uj). The lowest excitation in tower p cor- 
responds to a set of densely packed quantum numbers Ij 
shifted by p. From the definition of the momenta kj , this 
is equivalent to imposing a twist of wave- vector Q + 2pn, 
so that Aq(ui) — Aq +2p7T (ui). So all we have to do is cal- 
culate A^(uj), where Q can take values inside and outside 
the first Brillouin zone. Now, A°^(uj) has peaks at ener- 
gies £o + ju c 2 r K / L, j = 1,2,..., where eo is the energy 
of the lowest peak in the tower. This energy is equal to 
Eq +1 (Q) -Eg =e f + wu c (l + Oq)/L, where 



a Q 




(10) 



So AP-(uj) can be written as iA%S(oj — £o~ ju c 2-K/L). 
i i 

From Eq. (|9j), one can show that 

(1 + «o)(2 + o:a) ...(?' + «n) ~ n 

A 3 = ^ 92- Ri Y. 91 A % + 0(1/ L) 

Q j\ v 



(j+l/2 + ag/2)°a 
T(aQ + 1) Q 



(11) 



and 



A Q = -(Lshim)- a ^f(Q) + O(l/L) (12) 

where the dependence on size and filling is taken care of. 
If we put all together, we get 



1 



1 



^ 2iru c r(ag + 1) y 2ttu c sin irn 



uj — ef 



f(Q) (13) 



The function f(Q) for Q <C L satisfies the recursion re- 
lation 

f(Q + tt) = f (Q - tt) T ^1 2it)2 ^Ql- (14) 

y ' v ; r(-Q/2^) 2 v ; 

where /(7r) = 1 and in the interval from = to n it can 
be approximated within 0.01% by 

ln/(Q) re -0.3047 + 0.3248 Q 2 /n 2 - 0.0201 Q 4 /n 4 

Finally, due to the divergence of C a ^(Q) and D a ^(Q) 
at Q = 7r/2, A lhb (uj) is dominated close to Ef by 
A-k/2{oS). This function is itself dominated by A^ 2 (ui). 
which diverges with the exponent = ~3/8 at £p- So, 
A lhb (lj) cx \u>— £f\~ 3 ^ s near ef- The whole proof can be 
reproduced for B(lj) which also diverges as \uj — £f\~ 3 ^ s - 
Let us note that this power law cannot be clearly iden- 
tified numerically with large but finite systems because 
there are other subdominant diverging contributions. 

This divergence of A lhb (uj) with exponent —3/8 is dif- 
ferent from the prediction of the Luttinger liquid theory 
according to which A lhb (lo) should vanish as \u> — £f\ u , 
where a, the exponent that also enters the momentum 
distribution function, is known to be equal to 1/8 in the 
U — > +oo limit. This apparent contradiction can be lifted 
as follows. First, we note that the exponent —3/8 can be 
explained very simply on the basis of the Green's function 
of the large U Hubbard model p[] which is given by 



G{x,t) 



^(kfrX — SFt) 



1 



(a: - u c i) 1/2 (x - u s <) 1/2 (x 2 - u 2 c t 2 )V 16 
If we set u s = 0, then A{uj) behaves like 

p i(ui— ep)t 



dtG(x = 0,t)e 11 



dt- 



t 5/S 



[uj - £ F y 



3/8 



while, if u s is finite, it behaves like 



dlGir = ()./)<-' - / rft e ^ 9/ £ g F)t - (UJ-£ F ) 1/8 



However, in both cases, rik behaves like 



3 



/r p i(k-k F )x 
dxG(x, t = 0)e tkx ~ dx—^j- 8 \k-k F \ 1/8 

So, as long as u s > 0, or equivalently U < +00, the 
behaviour very close to £f is given by (oj — Ef) 1 ^ , but 
this behaviour is limited to frequencies (u> — £_f) < u s ~ 
t 2 /U. In the limit U — ■*■ +00, this domain shrinks to a 
single point, ui — ep, and the spectral function diverges 
as (w — £i?)" 3 / 8 . If {7 is large but finite, there will be 
a peak in the spectral function at lo ~ t 2 /C7 reminiscent 
of this divergence. This seems to be consistent with the 
Monte Carlo results of Preuss et al. |l3| which show an 
increase of the spectral functions close to the Fermi level. 

The implications for the experimental observation of 
the Luttinger liquid behaviour are rather dramatic. To 
be able to see a difference from a step function in spite 
of the finite resolution, a should be large enough, which 
means one should consider strongly correlated systems. 
But the present calculation shows that the range of valid- 
ity of the asymptotic law (ui — £f) q becomes very small 
if the correlations are too big. Whether realistic models 
of one-dimensional conductors with intermediate values 
of the repulsion terms can lead to measurable effects re- 
mains to be seen. 
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FIG. 1. (N+l)C a , N (Q) (open symbols) and (N-l)D a , N {Q) 
(solid symbols) for cluster sizes N = 10, 14, 18, 22, 26 



FIG. 2. A lhb (uj) and B(u>) for the quarter-filled U -> +00 
Hubbard model with L=260 and N=130. 
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